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' We investigate hard-thermal-loop (HTL) corrections to the final lepton asymmetry in lep- 

togenesis. To this end we solve the Boltzmann equations with HTL-corrected rates and CP 
asymmetries, which we calculated in paper I of this series. We pay special attention to the 
influence of the two leptonic quasiparticles that arise at non-zero temperature. We include only 
decays and inverse decays and allow for the lepton modes to be either decoupled from each 
other, or to be in chemical equilibrium by some strong interaction, simulating the interaction 
with gauge bosons. In two additional cases, we approximate the full HTL lepton propaga- 
tors with zero-temperature propagators, where we replace the zero-temperature mass by the 
(— I ' thermal mass of the leptons (T) or the asymptotic mass V2me{T). We compare the final 

lepton asymmetries of the four thermal cases and the zero-temperature case for zero, thermal 
and dominant initial neutrino abundance. The final lepton asymmetries of the thermal cases 
differ considerably from the vacuum case and from each other in the weak washout regime for 
\jr^ ' zero initial neutrino abundance and in the intermediate regime for dominant initial neutrino 

CO , abundance. In the strong washout regime, the final lepton asymmetry can be enhanced by a 

' factor of two in the case of strongly interacting lepton modes. 



1 Introduction 

The question of the origin of all things was always essential to mankind and has driven them to 
search for answers in science, among others. Physics as the science of nature and within physics, 
^ . cosmology as the science of the order and the evolution of the universe, address this question and 
^ I have their own formulation of it. What is the origin of the matter that is the building block of all 
things we observe, including ourselves? 

The matter in nature consists of leptons and the much heavier baryons, which are in turn made 
up of quarks. According to the standard model of particle physics (SM), matter particles, quarks 
or leptons, can only be created in pairs together with their antiparticles, that is, antiquarks and 
antileptons, at least in perturbation theory. If we assume that the early universe was indeed without 
form and void , that is in the language of particle physics, there was no excess of one particle species 
over the other, there would have to be an equal amount of particles and antiparticles today. More 
specifically, since annihilation of particles and antiparticles proceeds at fast rates, no structures 
like atoms, molecules, galaxies, stars, planets, DNA, cells and finally living organisms could have 
formed and we would observe a universe populated almost exclusively by photons and the slowly 
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interacting neutrinos. This scenario is obviously not realised. If we believe in inflation, we cannot 
assume a sizeable matter-antimatter asymmetry as an initial condition of the universe because that 
asymmetry would be diluted by inflation, not even mentioning the highly unsatisfactory character 
of such an approach from a scientific point of view. Therefore we have to employ a baryogenesis 
theory, a mechanism that creates a baryon asymmetry dynamically and explains the value of 



This value has been inferred from the WMAP seven-year cosmic microwave background (CMB) 
anisotropy data where ub, ti§, and n-y are the number densities of baryons, antibaryons, and 
photons, respectively, and the subscript implies present cosmic time. The value agrees with the 
abundance of light elements inferred from big bang nucleosynthesis. 

Leptogenesis [2] is a very attractive baryogenesis theory, since it simultaneously explains the 
creation of the baryon asymmetry and the smallness of neutrino masses via the seesaw mecha- 
nism [SHS]. We add three heavy right-handed neutrinos Ni to the SM, which are assumed to have 
rather large Majorana masses Mi, close to the scale of some possibly underlying grand unified the- 
ory (GUT) [9], -Egut ~ lo^^---^^ GeV. The interaction with the SM neutrinos suppresses their mass 
when we integrate out the heavy neutrinos. In the early universe, the heavy neutrinos decay into 
leptons and Higgs bosons and create a lepton asymmetry, which is lateron converted to a baryon 
asymmetry by the anomalous sphaleron processes (TOlttl] • The three Sakharov conditions [12] that 
are necessary for a baryogenesis theory are fulfilled, that is lepton number L and B — L are violated, 
CP symmetry is violated in the decays and the decays can be out of equilibrium. 

Ever since the development of the theory 25 years ago, the calculations of leptogenesis dynam- 
ics have become more refined and many effects and scenarios that have initially been neglected 
have been considerecH. Notably the question how the hot and dense medium of SM particles in- 
fluences leptogenesis dynamics has received increasing attention over the last years |14H21j . At 
high temperature, particles show a different behaviour than in vacuum due to their interaction 
with the medium: they acquire thermal masses, modified dispersion relations and modified helicity 
properties. All these properties can be summed up by viewing the particles as thermal quasipar- 
ticles with different behaviour than their zero-temperature counterparts, much like the large zoo 
of single-particle and collective excitations that are known in high density situations in solid-state 
physics. At high temperature, notably fermions can in the hard-thermal-loop-limit (HTL) occur in 
two distinct states with a positive or negative ratio of helicity over chirality and different disper- 
sion relations than at zero temperature, where these dispersion relations do not break the chiral 
symmetry as a zero-temperature mass does. 

Thermal effects have been considered by references |14H21j . Notably reference [15] performs 
an extensive analysis of the effects of thermal masses that arise by resumming propagators using 
the HTL resummation within thermal field theory (TFT). However, the authors approximated the 
two fermionic helicity modes with one simplified mode that behaves like a vacuum particle with its 
zero-temperature mass replaced by a thermal mas^. Due to their chiral nature, there are serious 
consequences to assigning a chirality breaking mass to fermions, hence the effects of abandoning this 
property should be examined. Moreover, it seems questionable to completely neglect the negative- 
helicity fermionic state which, according to TFT, will be populated at high temperature. We 

^For an excellent review of the development in this field, we refer to reference [13| . 

^Moreover, an incorrect thermal factor for the CP-asymmetry was obtained, as has been pointed out in refer- 
ence 122]. 
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argue in this study that one should include the effect of the fermionic quasiparticles in leptogenesis 
calculations and possibly in other early universe dynamics, since they behave differently from zero- 
temperature states with thermal masses, both conceptually and regarding their numerical influence 
on the final lepton asymmetry. We do this by analysing the dynamics of a leptogenesis toy model 
that includes only decays and inverse decays of neutrinos and Higgs bosons, but takes into account 
all HTL corrections to the leptons and Higgs bosons, paying special attention to the two fermionic 
quasiparticles. In a slightly different scenario, we assume chemical equilibrium among the two 
leptonic modes, thereby simulating a scenario where the modes interact very fast. As a comparison, 
we calculate the dynamics for two models where we approximate the lepton modes with ordinary 
zero-temperature states and modified masses, the thermal mass mi{T) and the asymptotic mass of 
the positive- helicity mode, \/2m^{T). 

This paper is the second part of a two-paper series, where we have calculated the HTL- 
corrections to CP-asymmetries in the first part ^23j . The topic of this work is solving the Boltzmann 
equations with HTL-corrected rates and CP-asymmetries. It is structured as follows: In section [2l 
we briefiy review the imaginary time formalism of thermal field theory (TFT) and discuss the hard 
thermal loop (HTL) resummation. In section [3l we summarise and present our previous calcu- 
lations of interaction rates and CP-asymmetries in references [21] , [25J and [23] . Section [J] deals 
with the evaluation of the Boltzmann equations. We derive the equations and compare our four 
thermal scenarios, wich are the decoupled and strongly coupled two- mode approach and the one- 
mode approach with thermal and asymptotic mass, to the vacuum case. We show the evolution 
of the abundances for three different initial conditions for the neutrinos, that is, zero, thermal and 
dominant abundance. We explain the dynamics of the different cases in detail and find consid- 
erable differences both of the thermal approaches to the vacuum case and of the two-mode cases 
to the one-mode cases. We summarise the main insights of this work in the conclusions and give 
an outlook on future work and prospects. In appendix |Bl we derive Boltzmann equations at zero 
temperature, while in appendix [Cl we explicitly perform the subtraction of on-shell propagators 
for our cases. 

2 Propagators at Finite Temperature 

When going to finite temperature [26] , one has to employ ensemble weighted expectation values of 
operators rather than the vacuum expectation values, so for an operator A we get 



There are two formalisms for calculating Green's functions at finite temperature, the imaginary 
time formalism and the real time formalism. Both are equivalent and we employ the imaginary 
time formalism, where the /cQ-integration is replaced by a sum over discrete energies, the so-called 
Matsubara frequencies. 

Naive perturbation theory at finite temperature can lead to serious conceptual problems, such 
as infrared divergent [271[28] and gauge dependent j29|[30] results and results that are not complete 
to leading order. In order to cure these shortcomings, the hard thermal loop (HTL) resummation 
technique has been invented [3HI32j. One distinguishes between hard momenta of order T and 
soft momenta of order gT, where g is the coupling constant of the corresponding theory. In a 
strict sense, this is only possible in the weak coupling limit where 5 <C 1. If all external momenta 
are soft, then the bare thermal propagators have to be replaced by resummed propagators. The 
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Figure 1: The resummed scalar propagator. 

self-energies that are resummed are the HTL self-energies, for which all internal momenta are hard. 
For a scalar field with a HTL-self-energy 11, the resummed effective HTL-propagator A* follows 
from the Dyson-Schwinger equation in figure [1] as 



iA* =iA + iA(-in)iA + ... 
i i 



A-i - n K^-ml-U' 



(3) 



where A is the bare propagator, K the momentum and mg the zero-temperature mass of the scalar. 
The dispersion relation for this effective excitation is given by the pole of the propagator as 

k'^ = e +ml + U, (4) 

so we get an effective mass of m^g = mg + m|, where the thermal mass of the scalar is given 
by the self-energy, which is proportional to gT, m| = 11 oc {gT)"^ . It is possible to neglect the 
zero-temperature mass if ms ^ rriQ. 

For fermions with negligible zero-temperature mass, the general expression for the self-energy 
in the rest frame of the thermal bath is given by [33] 

nP) = -a{P)^-b{P)i, (5) 

where = (1,0,0,0) is the four- velocity of the heat bath. The factors a and h are given by 

^(^) = i^[^'tr(7o5])-potr(/'S)] . 

In the HTL limit, the traces are given by [26] 

Ti = tr =4m|, 

^2 = tr (70S) = 2m^-ln — , (7) 

p Po-p + ie 

where the effective thermal fermion mass mp oc gT depends on the interaction that gives rise to 
the fermion self-energy. 

The resummed fermion propagator is then written as 

S*{K) = i — — . (8) 

^ ^ ^-SHTL(i^) ^ ' 
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Figure 2: The two dispersion laws for fermionic excitations compared to the standard dispersion 
4- 



relation 0^ = /c^ + mP, 



It is convenient to rewrite this propagator in the helicity-eigenstate representation |34y35j. 

S*iK) = \a+{K){jo - k • 7) + \a.{K){jo + k • 7), (9) 



where k = k/fc, and 



(10) 



This propagator has two poles, the zeros of the two denominators A±. The poles can be seen 
as the dispersion relations of single-particle excitations of the fermions that interact with the hot 
plasma, 

ko = 0±ik). (11) 

We have presented an analytical expression for the two dispersion relations making use of the 
Lambert W function in reference |24) . The dispersion relations are shown in figure [2l 

Note that even though the dispersion relations resemble the behaviour of massive particles and 
= mp for zero momentum k, the propagator S*{K) ([9]) does not break chiral invariance like a 
conventional mass term. Both the self energy Ti{K) ([5]) and the propagator S*{K) anticommute 
with 75. The Dirac spinors that are associated with the pole at ko = 0+ are eigenstates of the 
operator (70 — k • 7) and they have a positive ratio of helicity over chirality, = + 1. The spinors 
associated with ko = 0~, on the other hand, are eigenstates of (70 + k • 7) and have a negative 
helicity-over-chirality ratio, = — 1. At zero temperature, fermions have = + 1. The introduction 
of a thermal bath gives rise to fermionic modes which have = — 1. These modes have been called 
plasminos since they are new fermionic excitations of the plasma and have first been noted in 
references I33ll36l. 
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We can introduce a spectral representation for the two parts of the fermion propagator pO|) [37] , 

AJK)=r,.^i^. (12) 

where the spectral density p±{0,k) [MIEH] has two contributions, one from the poles, 

p^^^{0,k) = Z±{0,k) {0 - 0±{k)) + Z^{0,k) {0 + 0^ik)) , (13) 
and one discontinuous part, 

pf'i^^ = TT, — ^ — . ''"t;::2 — — ^ >< ^(^' - ^ (14) 



im|(A;=F0) 

{k{0^k)-m\\Q<i{x)^Qi{x)\]^ ^ \^Txm\{\^. 



X] 



where x = 0/k, 9{x) is the heaviside function and Qq and Qi are Legendre functions of the second 
kind, 

1 X + 1 

<3o(2;) = -ln -, Qi{x) = xQo{x) - I. (15) 

2 j; — 1 



The residues of the quasi-particle poles are given by 

0l{k) - e 

2 m 



Z±{0, k) = ^'^ ' ^^^6^6 Z+ + Z_ = 1 . (16) 



"F 

One can describe the non-standard dispersion relations 0± by momentum-dependent effective 
masses m±{k) which are given by 

m±{k) = ^0'i{k) - A:2 = v^2 Z(0, k) nip ■ (17) 

These masses are shown in figure [3l 

Considering gauge theories, one might also have to use HTL-corrected effective vertices that 
are related to the propagators by Ward identities [26]. We do not consider these vertices since we 
are only looking at Yukawa vertices. In the HTL framework, it is sufficient to use bare propagators 
if at least one of the external legs is hard. However, it is always possible to resum self-energies 
and thus capture effects which arise from higher-order loop diagrams and take into account the 
appearance of thermal masses and modified dispersion relations in a medium. In fact, since the 
effective masses we encounter do typically not satisfy the condition rngg <C T but are rather in the 
range rUes/T ~ 0.1 - 1, the effect of resummed propagators is noticeable even when some or all 
external momenta are hard. In summary, we always resum the propagators of particles that are 
in equilibrium with the thermal bath, which are the Higgs bosons and the leptons in our case, in 
order to capture the effects of thermal masses, modified dispersion relation and modified helicity 
structures. This approach is justified a posteriori by the sizeable corrections it reveals, similar to 
the treatment of meson correlation fuctions in reference |39] . 

In leptogenesis, the leptons and Higgs bosons acquire thermal masses that have been calculated 
in references [SSlIMllinilll] and are given by 
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Figure 3: The momentum-dependent effective masses m±. 



The couphngs denote the SU(2) coupUng (72, the U(l) couphng gy, the top Yukawa couphng yt and 
the Higgs self-coupling A, where we assume a Higgs mass of about 115 GeV. The other Yukawa 
couplings can be neglected since they are much smaller than unity and the remaining couplings are 
renormalised at the first Matsubara mode, 2ttT, as explained in reference [15] and in reference |42j 
in more detail. The heavy neutrinos A^i do acquire a thermal mass, but since the Yukawa couplings 
are much smaller than unity, this effective mass can be neglected compared to the zero-temperature 



mass. 



3 HTL Corrections to Decays and CP- Asymmetries 
3.1 Decay and inverse decay rates 

The additional terms for the right-handed neutrinos in the Lagrangian are 

C = i N.d^j^'Ni - liaNiireabt) -^Yl MiNiN^ + h.c. , (19) 

i 

where the Higgs doublet (j) is normalised such that its vacuum expectation value (vev) in 

i^) = (°) (20) 

is f ~ 174 GeV and lia is the Yukawa coupling connecting the Higgs doublet, the lepton doublet 
and the heavy neutrino singlet. The indices a and b denote doublet indices and e^b is the two- 
dimensional total antisymmetric tensor that ensures antisymmetric SU (2)-contraction. 

We have discussed the HTL corrections to neutrino decays Ni — )• HL in detail in reference |24j . 
When the temperature is so high that > Mi, the neutrino decay is kinematically forbidden 
in the HTL-approximatior|l, but the decay of Higgs bosons into neutrinos and leptons becomes 

^It has been shown in reference [17] , that the decay is still allowed if one considers the effect of collinear external 
momenta. 
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Figure 4: A^i decay via the optical theorem with dressed propagators denoted by a blob. 



possibl^ [23]. We have briefly shown the Higgs boson decay rate in reference [,25j and explain the 
calculation in more detail in a separate work |23) . The basic idea is to calculate the rate via the 
optical theorem by cutting the neutrino self-energy with resummed propagators as in figure |H Using 
this method, we describe the particles that are affected by thermal corrections, the Higgs boson 
and the lepton, by thermal propagators, whereas the external particle is not affected thermally 
since the couplings are small. According to finite-temperature cutting rules [431144] . the interaction 
rate for a neutrino with momentum P reads 

r(P) = --^ trii^ + Ml) Im S(po + ie, p)]- (21) 

Integrating over neutrino momenta p, we can write the decay density in a familiar form that is 
corrected for the statistical distribution of the particles, 

7(iVi^i7L) = J dpM,dpHdpU27r)^5\pN,-pL-PH)\MH\^f^\{l-fl''){l + fj,^) (22) 

and in the same way for the Higgs boson decays 

^{H^N^L) = J dpN,dpHdpU27T)^5^{pN,+PL-PH)\MH\'{l-fNV('^-fr)fH\ (23) 

where dp = d'^p/[(27r)^2£'], we have Ejj = p'j^ + ?n.^, = ^h{PL) and /i = ±1 denotes the helicity- 
over-chirality ratio of the lepton. The inverse processes HL — ?■ Ni and NiL ^ H can be written in 
the same way with the appropriate statistical factors. The matrix element for neutrino and Higgs 
boson decays turns out to be the same, like at zero temperature, 

\Mh{P.K)\' = A{in)iiP^K^ = A{in)iiZhUJh{PG - /ip • k), (24) 

where we have introduced a chirally invariant four-momentum = Zh[h)0h{k){l,h'k) for the 
lepton and k = k/Zc. From this matrix element, one can derive a multiplication rule for the HTL 
lepton spinors. 



and the antiparticle spinors. 



iK)uf (K) = Z±u:± T k • 7) , (25) 



vf{K)vf{K) = -Z±^±(7o ± k • 7) . (26) 



The lepton decay is not possible, since > mi for all temperatures 
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Figure 5: The decay densities for the neutrino and the Higgs boson decay. We show the one-mode 
approach with the thermal mass as and with the asymptotic mass as 7^^^; Also the T = 
rate 70 and our two modes 7-1- . The temperature thresholds are explained in the text. 



These rules make it easy to calculate processes that involve HTL-corrected leptons as external 
particles. 

In figure m we compare our consistent HTL calculation to the one- mode approximation adopted 
by reference [15], while we add quantum-statistical distribution functions to their calculation, which 
equals the approach of using an approximated lepton propagator l/(-?t — mi) [45]. In addition, 
we show the one- mode approach for the asymptotic mass \plmi. We evaluate the decay rates for 
M\ = 10^*^ GeV and normalise the rates by the effective neutrino mass m\ = (/^Z)iif^/Mi, where 
V = 174 GeV is the vacuum expectation value of the Higgs field. This effective mass is often taken 
as 777-1 = 0.06 eV, inspired by the mass scale of the atmospheric mass splitting. 

The decay densities are analysed in detail in reference [23] . Summarising, we can distinguish five 
different thresholds for the thermal decay rates. Going from low temperature to high temperature, 
these are given by the following conditions: 



rr,N 


Ml 


= \pi mn + 


rpN 


Ml 






Ml 


= , 


-'o 




= 1711 + Ml, 


rp4> 
^ + 




= V2mi + Mi 



These thresholds correspond to three different thermal lepton masses, the asymptotic mass of the 
(+)-mode, ^/2m£, the naive thermal mass 777^, which is the effective mass for zero lepton momentum. 
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Figure 6: The CP-asymmetry in neutrino decays. The graph in the middle is the self-energy 
contribution, the graph on the right the vertex contribution. 




Figure 7: The CP-asymmetry in Higgs boson decays at high temperature. Again, the graph in the 
middle is the self-energy contribution, the graph on the right the vertex contribution. 



and the vanishing asymptotic mass of the (— )-mode. 
3.2 CP-asymmetries 

The CP-asymmetry in neutrino decays at zero temperature is defined as 

° T{N (Pi) + r{N ^ 0£) ■ ^ ' 

At finite temperature, we have to calculate the CP-asymmetry via the integrated decay rates, 

where we define the CP-asymmetry for each lepton mode, denoted by h. The CP-asymmetry arises 
as an interference between tree level and one-loop diagrams in neutrino and, at high temperature, 
Higgs boson decays, shown in figures [6] and [71 In order to calculate the interference, one has to 
take the imaginary part of the relevant one-loop diagram. At zero temperature, this can be done 
by cutting through the diagram and determining the discontinuity via the optical theorem and the 
Cutkosky rules. Kinematically, it is only possible to put the lepton propagator and the Higgs boson 
propagator in the loop on-shell, so there is one possible cutting for the self-energy graph and one 
possible cutting for the vertex correction graph. 

At finite temperature, there exist cutting rules for the RTF. In the ITF, it is possible to isolate 
terms with certain momentum relations that correspond to certain cuttings, while a direct relation 
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to the RTF cutting rules is not straightforward. Since we can exchange energy with the heat bath, 
also terms that correspond to cutting through the in the loop are possible [15 11211 - 123] . In the 
hierarchical limit M2 3> Mi, which we assume here, these terms are suppressed, so again the cuts 
through the lepton and the Higgs boson in the loop survive. Due to the two possibilities for lepton 
and Higgs doublets in the loop, the CP-asymmetry from the self-energy graph is exactly twice as 
large as the CP-asymmetry from the vertex correction graph in the hierarchical limit. We can 
therefore distinguish four contributions to the CP-asymmetry, taking into account the two lepton 
modes in the loop and two modes for the external leptons. 

The differences in neutrino decay rates, that is A7 = 7(A^ — )• (j)tj — j{N — t- (j)£), are shown 
in figure [8] and compared to the two one-mode approaches with masses m£ and \pl'mi. The 
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Figure 8: The CP-asymmetries in neutrino decays normalised by the CP-asymmetry in vacuum 
and the total decay density in vacuum, A7/(7Q°*eo). We choose Mi = 10^^ GeV and M2 ^ Mi. 
The term A'j^j^h^ denotes the difference between the decay rate and its CP conjugated rate, which 
is proportional to the CP-asymmetry. Here, hi denotes the mode of the external lepton, while /12 

denotes the mode of the lepton in the loop. For example, A7_| = 7(A^ — )• 0^+) — — )• 0^+), 

where a minus- mode lepton is present in the loop. A^jm^ and ^7^^^ denote the rate differences 
for the one- mode approach with a thermal mass m£ and an asymptotic thermal mass \f2mi. 

same is shown for Higgs boson decays in figure [9l We discuss these CP-asymmetries in detail in 
reference 1231. 



11 



o 

IS 
o 



10000 








100 








1 

0.01 






y ^. 




/ 
1, 

// 
a 

■-f 

( 

1 




0.0001 

10-6 


,' 


^'V2mi, 

A7++ 

A7+-/-+ 

A7 


10-8 






. .1 . . 



10 



100 



T/M 



Figure 9: The CP-asymmetries in Higgs boson decays normalised by the CP-asymmetry in vacuum 
and the total decay density in vacuum, A7/(7Q°*'eo), where the asymmetries A7 are explained in 
figure El We choose Mi = 10^° GeV and M2 > Mi. 



4 Boltzmann Equations 

We calculate the Boltzmann equations for leptogenesis. We include decays and inverse decays 
involving neutrinos, leptons and Higgs bosons. We neglect all scatterings, expect the on-shell 
contribution of the AL = 2 scatterings, which we have to take into account for consistency rea- 
sons. We take into account thermal dispersion relations, but assume distributions close enough 
to equilibrium that we can use Boltzmann equations. For the distribution functions, we use the 
full quantum statistics, that is, Fermi-Dirac and Bose-Einstein statistics, but assume the kinetic 
equilibrium approximation /j = rii/n^^f^'^. It has been shown by reference [IB] that this is a good 
approximation. The Boltzmann equations at zero temperature are derived in appendix [Bl 



4.1 Low Temperature 
Neutrino evolution 

The formulation of the Boltzmann equations we need is given by equation (jl34p in appendix |Aj 
The equation for the evolution of the lightest right-handed neutrino reads 



+ 7(7Vi 



) + l{Ni 
+ liNi - 



) - 7(0^- 



> Ni) - ^{^i, 



Ni)] 



(30) 



where i± denote the two lepton modes. We neglect scatterings since they are of higher order in the 
coupling constant. We will from now on omit the subscript 1 for the neutrino and write A^. The 
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CP- asymmetry in the matrix element is not relevant for neutrino decay, so we calculate the matrix 
element, which is the same for the above processes and define 

\Mlf = \M{N ^ HL±)\^ = \M{HL± ^ N)\'^ , (31) 

where now HL± denotes the sum of leptons and Higgs doublets, i± and <j), and their charge 
conjugated states i± and (p. The matrix elements are, however, different for the different lepton 
modes and also the momentum-conserving delta functions differ from each other. The subscript it 
means (+) or (— ), not the sum. When summing an expression A± that is dependent on the kind 
of lepton dispersion relation over the lepton modes, we write ^_|_ A±. We have 

For each of the two lepton modes, we have now 

7(iV ^ HL±) - 7(FL± ^N) = j dpN,dpL±dpH{2T^f5\pN-PH-PL±) 

X \M{N ^ HL±)\^ fN{l + fH){l- fL±) 
- \M{HL± ^ Nf (1 - fN)fHfL±\ . (33) 
The term in square brackets in equation ([33]) reduces to 

\Ml\^ [/7V(1 + Mil - fL±) - (1 - fN)fHfL±] = \Mlf [CN^HL± " ChL±-,n] , (34) 

where 

CN^HL± = /Ar(l + /h)(1 - /l±) , 

CHL±^N = (1 - fN)fHfL± ■ (35) 

Throughout this section, we make the kinetic equilibrium assumption, that is, the phase space 
densities can be written as 

f^ = = ^^fr, (36) 

where Xi = nj/n^^, and n^'^ and f^"^ are the equilibrium number densities and distributions. For 
the neutrino evolution, we can assume that the Higgs bosons are in equilibrium since they couple 
very strongly to the thermal bath, fn = fj^- The lepton distributions of the two modes are, 
strictly speaking, out of equilibrium since leptons and antileptons are created asymmetrically. 
However, the leptons are much closer to equilibrium than the neutrinos, so the neutrino evolution 
is not influenced by the lepton asymmetry. Therefore we approximate the lepton densities with 
their equilibrium density, fL± = fi±- We will relax this assumption in the section on the lepton 
asymmetry evolution. 
Using the relation 

r^ii + /^^)(i - fii) = (1 - /^^)/r±/^^ (37) 

we can write 

CN^HL± = XnQO- + )(1 - /l±) = {XN - XNQ)fHfL± ' 

CHL±^N = (1 - XNf^)f^fl± , (38) 
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and 

CN^HL± - CN^HL± = {xn - (39) 

The decay densities are 

7(A^ ^ HL±) - -f{HL± ^ N) = j dpNdpL±dpHi2iT)'^6'^{pN - Ph - Pl±) 

= {xN - 1)7d±, (40) 

where 

7^± = I dpNdpL±dpH{27r)'S\pN-PH-PL±) |Xi|'/K±- (41) 

Note that 7^_|_ is not the same as the equihbrium decay density in equation ()22|) . but differs from 
the latter through the thermal factor JhIl- It is an effective decay density, which enters the 
Boltzmann equations. The Boltzmann equation for the neutrinos reads 

where 7^ = 7+ + 7^, or, in analogy to equation (|146p in appendix [Bl 

^ = -D^in^-n'^^), (43) 
dz 

where 



^^ = %T^ (44) 
rij^ Hz 

and we have used H{Mn) = Hz^. Most conveniently, the number densities are normalised by the 
entropy density s in order to factorise their dependence on the expansion of the universe. The 
entropy density scales as 

s = 9.^T^ (45) 
where (7* counts the total number of effectively massless degrees of freedom and is defined as 

E ».(f)*+^ E f,(f)\ (46) 

j=bosons j=fcrmions 

where i denotes species with mass rrii <^T and the factor 7/8 arises from the difference in Fermi 
and Bose statistics [17]. At the temperature of leptogenesis, all SM particles have negligible masses, 
so 5* = 106.75. We define all number densities in terms of the entropy density as 

= ^, (47) 
s 
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then 



Xi = ^. (48) 



The Boltzmann equation reads 

'^^'^ " (xN-lhE, (49) 



dz sHi 
where 



and the Planck mass is 

Mpi = 1.221 • 10^*^ GeV. (51) 

Lepton asymmetry evolution 

We set up evolution equations for the two different lepton modes separately and define the phase 
space density of the lepton asymmetry in the respective mode as 

fch = feh - fih- (52) 

where h = ±1 denotes the helicity-over-chirality ratio of the leptons. The final lepton asymmetry 
is then n^" = n^"^ + after evaluating the Boltzmann equations for each mode separately. The 
Boltzmann equations for leptons and antileptons read 

^^-^{7(4.0-^) -7(^-4.0) 

+ E b^^^l*^ ^ ^"^2 0) - ^^^^^^ ^ ^hA)] } , 

^ = -^{7(4.^-iV)-,(^^4.^) 

+ E [7(4i<A ^ ihj) - l{ih,(t> ^ 4i0)] }, (53) 

where we have {hi, /12) = ±1 to account for the second lepton involved in the scatterings and we 
have only included AL = 2 scatterings since the other two-by-two scatterings that involve only N, 
i and (j) are negligible |15j . For the evolution of the lepton asymmetry, we have 

= - 7^|7(4i0 ^N)- 7(4i0 -^N)- 7(iV ^ ih,<P) + j{N -> 4i0) 
dz Hz L 

+ E b^^h,^ ^ ^hj) - 7(4i<^ ^ 4^0) + 7(4^0 ^ 4i<^) - 7(42<^ ^ 4i</')] } • (54) 

At leading order in the couplings, the AL = 2 scatterings are computed at tree level and are 
consequently CP-conserving. However, from these scatterings we must subtract the CP-violating 
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Figure 10: The s-channel contribution to the AL = 2 scattering 



contribution where an on-shell A^i is exchanged in the s channel, shown in figure [TUJ This is 
because in the Boltzmann equations the process is already taken into account by inverse decays 
with successive decays, £0 — )• — )• [l5 t H8 l H9]. We must therefore replace the scattering rate by 
the subtracted rate. 



7(4i0 ^ th^<^) ^ 7'"'(4i</' ^ 420) = liihA ^ 420) - 7°"~''''''(4i0 ^ 4^0) (55) 

and 'y{ihi4' ~^ ^/i20) accordingly, where ^o'^-'^hsll jg ^^iq on-shell contribution. The Boltzmann 
equations then read 



dz 



+ E [7™'(4,0 ^ 420) - f ""(4,0 ^ 420) 

ft2 

+ 7^""(420 ^ 4i0) - 7'""(420 ^ 4i0)] } . (56) 



It is convenient to define a CP- asymmetry in neutrino decays on amplitude level as 

^ |.M(iv^ 04)1' -|.M(iv^ 04)1' 
"'^ |A^(Ar^04)|2+ |A^(A^^04)|' 



and \MiN (t)th)\^ + |X(Af ^ 04)1 = l-^oftP, we write 

1 -I- 

|A^(iV ^ 04)1' = |A^(04 ^ iV)|' = lA^ohl' ^ 



2" 



1 _ jv 

|A^(iV ^ 04)1' = |A^(04 ^ iV)|' = \Moh\' ■ (58) 

It is useful to write the decay rates for the above 1^2 processes as in section 14.11 

7(process) = j \{dp,{27:)H^ Pj) c(process), (59) 

where pj denotes the relevant momenta pN, Pih = Plh and p^ = p^ and 5^{Y^Pj) the momentum 
conservation 5^{pN —peh ~P(j>)- The information about the specific process is encoded in c(process) 
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and we have 

CN^,^4> = \M{N^ 04)1 Viv(l - hh){l + U) 
a^^^N = IMi^lh ^ A^)|' (1 - fN)fehU 
c^^l,-^ = \M{N^ fh)\^fN{l - hh){l + U) 

Ci^f^N = ^ iV)|' (1 - fN)fehU. (60) 

Since we are looking at the lepton asymmetry, the lepton distributions have to be out of equihbrium, 

fi/ih = ^e/ihfeh^ 

fch — ^C.hf^hj 

hh + flH « Vtl (61) 

while the Higgs bosons can be assumed to be in equilibrium. 

As explained in appendix O the scattering rates can be written as 

= J dpNdpe,Ap42nf5HpN - pen. - P4>)€ 1-^° Qf^i^ " Q)- (62) 



so we 



defin^ 



ct' = 2e^\Ml\'f2C{l-Q) (63) 
and calculate the integrand for the right-hand side of the Boltzmann equation (|154p , 

c{N ^ ^h(t^) - cilhcP ^ N)-c{N ^ ih4>) + c(40 ^N) + = 

= xcnmr;" + xn^ - 2e^Qr^ {xm - 1) (1 - 2/^^) . (64) 

We can easily check that this term vanishes when the neutrinos are in equilibrium, xat = 1, and 
there is no previous lepton asymmetry, xch = 0. 
The Boltzmann equation reads now 

= -Jj^ [-^-yhleh [XN - 1) + — {-fwh + XNlNh)\ , (65) 

where = is defined in equation (j^T]) and 

lS^ = f dpNdp,Hdp^{27i)H\pN -P4.- Pm) \Mo\' f^'Qil - 2Q) 
iNh = [ dpNdpihdp^i27r)'^5'^{pN -P4>- Peh) I-MqI^ frnfi^' 



— / dpNdpehdp^{27r)'6\pN -P^- Peh)e^ l>^o| V;V,T(1 - 2/^^). (66) 



^Note that our factor c°"'' differs from reference 3fi]i where they have the out-of-equihbrium distribution (1 — /jv) 
instead of (1 — fj^)- However, as derived in appendix[C] we must employ f^, even if we had only one lepton mode, 
which also results in a Boltzmann equation for {£ — £) which is slightly different from reference [46] . 
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We see that the rates and the CP-asymmetries that enter the Boltzmann equations have slightly 
different thermal factors than the equilibrium rate in equations (j22p . which employs the factor 
fN(.l-mi + U) for iV decays. 
We may also write 



dYch 



dz sHi 
or, corresponding to equation (jl54p in appendix [Bj 



e^HlS. i^N - 1) + ^ (7^. + x^j^,)] (67) 



where 

Di 



4.2 High temperature 

As discussed in section 13.11 the neutrino processes N ■H' are forbidden when the thermal masses 
of the Higgs bosons and leptons become too large, that is, when > Mj\f. However, new processes 
with the Higgs as single initial or final state are then allowed, <j) o A^^. These are the dominant 
contributions to the neutrino and lepton evolution and they can be CP-violating as well, so they 
contribute to generating a lepton asymmetry. We derive the Boltzmann equations for this high 
temperature regime in the following. 

Neutrino evolution 

We derive the Boltzmann equation analogously to section W7\\ 

h 

We have 

j{NLh ^H)- j{H ^ NLh) = j dpNdpLhdpH{2^f5\pH - Pn - PLh) 

X \M{NLh^H)\^fNfLh(.l + fH) 

-\M{H^ NL^,)\' (1 - /^)(1 - fLh)fH\ ■ (71) 

The tree-level matrix elements are the same at high temperature for the Higgs-processes, 

just the kinematics differ. So we have 

\Ml\^ = \MiNLf, ^ H)\' = \M{H ^ NL^)\^ . (72) 
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Again, we assume the Higgs bosons and leptons to be in equilibrium. We write 

\Mo\^ [fNhhil + Jh) - (1 - /iv)(l - fLh)fH] = \Mo\^ [c{NLh ^H)- c{H ^ NL^)] 
Using the relation 

+ fi^) = (1 - - flDfZ 

we get 

c{NLh ^H)- c{H ^ NLh) = {xm - 1)(1 - fZ)fH- 
The Boltzmann equation then reads 

— (xjv - 1)7d> 



dz Hz 

AYn __ z 
dz sHi 

or 



{xn - 1)7£), 



dz 



= -D'^{nN - n^), 



where 7^ = 

ith = I dpNdpLHdpH{2^)'S'{pH -PN- PLh) \MI\' f^{l - fH) 



and 

rQ Hz 

Lepton asymmetry evolution 

The Boltzmann equations for leptons and antileptons read 
dneh, z 



+ J2 bi^hA ^ 4,-^) - 7(4. <^ ^ 4i0)] }, 



/12 

combined we get 
dnchi 



dz 



= - -^{^(^h.N ^^)- -i{lh,N ^<P)- 7((^ ^ th,N) + 7((/) ^ 4iA^) 



+ 

/i2 



[i{h,4> ^ 42 0) - 7(4i^ ^ 42 0) + 7(42 ^ 4i0) - 7(42 ^ 4i0)] } • 
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Figure 11: The li-channel contribution to the AL = 2 scattering — )• l^. 

At high temperature, there can be no on-shell neutrino in the s-channel of the AL = 2 scatterings, 
but there can be an on-shell neutrino exchange in the n-channel as shown in figure [H] Again, we 
need to subtract the AL = 2 rates since the n-channel on-shell neutrino exchange corresponds to 
a Higgs decay followed by an inverse decay, — )• INl — )• ^l. We replace 

7(4i0 ^ 4.<A) ^ 7"^^4,</' ^ Ik,4>) = 7(4, <A ^ 4.<^) - 7r'''^'^4,<A ^ 4.<^) (84) 



and get 
dn£fe, _ 
dz 



-i-{7(4iiV ^ 0) - 7(4i ^ </>) - 7('A ^ 4i A) + 7(0 ^ 4i A) 



+ 



^ [7"^^4i</> ^ 4,0) - 7'"^4i0 ^ 4,0) + 7^^4,0 ^ IhA) - 7'"^42<^ ^ 4i 

hi 

We define a CP-asymmetry in Higgs decays on amplitude level as 



(85) 



thus 





_ \M{4)- 


^ A4)|'- 


\M{4> - 


^ A4)|' 


(86) 




~ \M[4>- 


^ A4)|' + 


\M{4>- 


^ A4)r 


\MQ)- 


> A4)|' = 


|A^(4A- 


>0)f = 






\M{(\)- 


> A4)|' = 


|A^(4A^- 


^0)l' = 


X-e'*' 


(87) 



As explained in appendix [Cl the scattering rates are written as 



[7^^^4> ^ 4,0) - f "^40 ^ 4, 



E [7^"'(4,</' ^ 4,<^) - 7™'(4,<^ ^ 4,c 

y dpArdp£ft.dp</,(27r)'^(5^(pAr - - P'f>)^h I 



KJ'/2(l+/^)/^^ 
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so 

cr' = 26tK|'/,T(l + /r)/^^ (89) 
and we get for the right-hand side of equation ([85]) . 

c{lhN ^ (^)-c(4iV ^ 0) - c{4> ^ ihN) + c{cl) ^ £hN) + 4"^ 

= ^cnmr;" + ^^/^^) + 2et(l - /.T)/;*^ (x^ - 1) (1 - 2/^^) . (90) 
The Boltzmann equation reads now 



dnch 1_ 

dz ~ Hz 



A^tu i^N - 1) + ^ (7t. + ^ivTt)] > (91) 



where 



%h = J dp,vdmdp42^)V(p^-p<^+p,^)|A^^|V7(l-/,7)(l-2/^^ 
Tat/i = / dpArdp£/idp0(27r)^5'^(pAr -p^+pih) /a/^"^' 



4 / dpNdp,hdp^{2^Y5\pN -P^+P,h)ei \Ml\^f^{l - - (92) 

7 !, J 



Teh 

Analogous to equation (jl54p . we may also write 



^ = -4^i(-^ - n;^) - (< + Wt„x^)ncH, (93) 



where 







< = 


1 7^. 
/fz 2<^ 






2n^^ 



(94) 



Normalised by the entropy density, the equation reads 
dYch z 



-e. 



'hith (^^ - 1) + ^ (iwh + ^^7^) ] • (95) 



dz si/i L T 
4.3 Interacting Modes 

The Boltzmann equations in the previous sections were derived under the assumption that the only 
relevant interactions in which the leptons take part are the Yukawa interactions with Higgs bosons 
and heavy neutrinos, which have very small coupling constants A, while it is implicitly assumed that 
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gauge interactions keep the leptons and the Higgs bosons in equiUbrium. This scenario would imply 
that the two modes only interact with each other via intermediate neutrinos or Higgs bosons, where 
the distributions and also the asymmetries in each mode are to first approximation decoupled. In 
a more realistic model, the Icpton modes will couple to each other via the SU{2) and U{1) gauge 
bosons and in processes like i± — >■ Izs^B. While it is conceptually interesting to consider 
the case that the two modes are completely decoupled, it might be more realistic to study the 
scenario where the interactions between the lepton modes are fast enough to keep them in chemical 
equilibrium. 

Chemical equilibrium implies that for species that interact via processes a + h — )■ i + j, the 
corresponding chemical potentials are related as 

l^a + l^b = + l^j ■ (96) 

When the processes which create or annihilate the particles and antiparticles of some species arc 
fast, for example via a + a — ^ i + j, where i and j are in equilibrium and their chemical potentials 
vanish, then the chemical potentials of a and a behave as 

IJ-a + IJ-a = IJ'i + IJ'j = , 

=^ IJ-a = -fJ-a ■ (97) 

In order to derive the corresponding Boltzmann equation, we introduce a chemical potential /x/j 
for the lepton mode £h- For simplicity, we approximate the distribution with Maxwell-Boltzmann 
statistics, an approximation which is sufficient to derive the final Boltzmann equations. The dis- 
tribution functions in kinetic equilibrium arc 

fihik) = e-'^("'^-^'') , 

hh{k) - Mk) = e-^^^{e^^^ - e-^^'^) ^ , (98) 

for <C T. We assume chemical equilibrium between the plus- and the minus-mode, 

H+ = li- = Hi. (99) 

Moreover, we can make the approximation that the equilibrium densities are about the same since 
the thermal 0.2 T is too small to affect the momentum integration considerably in 

eq [ ^ k \ 

"^'^ " J (2^^^'^ ^ ^ ' 

nf^ « « n|^o , (100) 
where is the distribution for massless leptons. With these approximations, we have 

XC± = — eq- , (101) 



n 
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where the subscript it indicates that we sum over the two modes, contrary to its use in the previous 
sections. We can now add the Boltzmann equations for the two modes in equations ([67ll and (f95]) 
and arrive at 



dz sHi 

where 



A7± (xN - 1) + — ;— {'yw± + XN7N±) , (102) 



Yc± = Yc+ + Yc- , 
A7± = e^+7e+ + 67- 7e- , 
lw± = lw+ + iw- , 

1N± = 1N+ + IN- ■ (103) 

The factor 1/4 comes from the fact that xch = xc±/'^- Depending on the temperature regime, we 
either have to employ the Higgs boson or the neutrino rates in the Boltzmann equations. 



4.4 One-Mode Approximation 



As we did in section [3] for the decay rates and the CP asymmetries, we also employ the one- mode 
approach for the Boltzmann equations. The equations are derived in analogy to sections 14. H and 
and read 



dz 
dYc 
dz 



z 



sHi 

z 



{xN - 'i-hDme, 
^Ime (xN - 



1) + — {iWmi, + XNlNme] 



(104) 



where jDmi-, IWme, iNmi and A7m^ are the same as the rates defined in equations (jlT]) . (j66|) . ([79]) 
and (192p and one has to make the appropriate replacements for the matrix elements and the lepton 
dispersion relations of the one- mode approach for and \/2m^. 



4.5 Evaluation of the Boltzmann Equations 

We solve the Boltzmann equations for five different scenarios: 

1. the zero temperature case with Maxwell-Boltzmann statistics, 

2. the two-lepton-mode approach where the two modes do not interact with each other, 

3. the two-mode approach where the modes couple strongly to each other, 

4. the one-mode approach for a thermal mass m^, 

5. and the one-mode approach for an asymptotic thermal mass \/2m^. 

In the decoupled case, the lepton asymmetries for the plus- and the minus-mode evolve separately 
from each other. When solving the equations, one has to specify the initial conditions for the 
neutrino abundance and the lepton asymmetry. We assume a vanishing initial lepton asymmetry 
and distinguish between three cases for the neutrino abundances: 
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1. Zero initial abundance: this is the case, for example, when an inflaton field decays mostly 
into SM particles and not into the heavy neutrinos. 

2. Thermal initial abundance: this can be realised when some additional interactions keep the 
neutrinos in equilibrium at T Mi, for example via a heavy Z' boson related to SO(VS) 
unification [501. 



3. Dominant initial abundance: this is the case, for example, when an inflaton decays predomi- 
nantly into A^i. 

The coupling (A^A)ii, which enters the neutrino decay rate, is parameterised by the so-called decay 
parameter defined as 



where 



= (105) 

mi = (106) 

is the conveniently defined effective neutrino mass, which is of the order of the light neutrino mass 
scale, and 



Ml 



-l.lxlO^^eV (107) 

r=Afi 

is called equilibrium neutrino mass. In the language of these masses, the out-of-equilibrium condi- 
tion, T < H, corresponds to > 1 and is called strong washout regime. The case < 1 is called 
weak washout regime. 

We want to analyse the evolution of the neutrino abundance and lepton asymmetries for the 
weak and strong washout regimes and different initial abundances. To this end, we write the 
Boltzmann equations for the different scenarios in the form of equations (I43p , (I68p , (I78p and (I93p , 



-D{Ym - Y^ 

dYr 

-^ = eoD,{YN-Y^'')-{W + WNXN)Yc, (108) 



where 



Hi sY^^ ' H, eosY^"^ ' 

^ = ^2^' ^^ = ^2^1^- 

One usually refers to D^iYf^ — Y^) as source term since this term is responsible for the production 
of a lepton asymmetry. The term {W + Wnxm)Y^ is called washout term since it usually has the 
opposite sign as the source term and reduces the production of the lepton asymmetry. The terms 
D, D^, W and Wn are different for the different scenarios. Note that for the finite temperature 
cases, is not the same as D and there is an additional washout term Wn due to the quantum 
statistics. Our analysis closely follows the arguments and explanations in reference [l9] and the 
interested reader will find a comprehensive explanation of leptogenesis dynamics in the vacuum 
case therein. 
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Weak washout for zero initial abundance 



Let us start with the weak washout regime and zero initial abundance. We define a value Zeq by 
the condition 

YN{Zeq) = Y^\Z,^) . (110) 

For z <C 1, the neutrino abundance is negligible compared to Y^, 

(111) 

where is approximately constant for z <^ 1. The entropy density s is proportional to z~^ 
and is proportional to z~'^ in vacuum and for the Higgs boson decays at high temperature 
in the finite temperature cases. Thus D ~ in the vacuum case and D ~ const, at finite 
temperature. Neglecting y^iti^.^ and z™^^^'^^, the integration yields Yj\[ ~ zD{z)/3 ~ z^ for the 
vacuum case and Yat — zD ~ z for the finite temperature cases. We show the numerical results 
for K = 0.005 and zero initial abundance in figures [T2l and [T3l where these power laws for Yn{z) 
can be observed for z < 0.1. We also see that > y/^'"^ > Y^ ^ Y^, which reflects 

> 7^™^ > 7^ 7o- Between the thresholds and z^ , the finite-temperature abundances 
do not evolve much, which reflects that the decay rates are very low or vanishing in this regime. 
The neutrino abundance for the asymptotic mass yplrrn does not rise at all at high temperature, 
while Y^^ rises slightly between z^^ and ztj.^ ■, where the rate is non-zero. The two-mode rate 
7-1- is, though very suppressed, present over the whole threshold range between zt^ and z^ due to 
the minus-modes, so Y^ rises slightly. At low temperature, z > Zq^, the neutrino abundances of 
the different scenarios are very close to each other, since for z > 2, the rates are very close to the 
vacuum rate, 75'^jv — '^I'^'^'^/^o — 7o- In this regime, Y^ is much larger than since the 
coupling is too small to keep the abundance close to equilibrium. 

Having outlined the evolution of the neutrino abundance at high temperature, we can un- 
derstand the evolution of the lepton asymmetry at high temperature. The washout term in the 
Boltzmann equations proportional to Yc is much smaller than the source term D^(Y]sf — Y^), since 
Yc/^o ^ Y^, and we have for z <^ z^q, 



Yl=o^-zeo^. (112) 



In the high temperature regime, the lepton asymmetry is negative and follows the neutrino abun- 
dance in its absolute value, 

Yc{z) -eo^Y^iz) = -^Ym{z) . (113) 

For the vacuum case, = D and Yc/cq ~ —Y^, while for the finite temperature cases, 
D^/D = A7/(eo7D) ~ 0(10^)11 Both behaviours can be observed in figures [T2] and [T3l The 



*We saw in section IT^ that A7/(eo7D) ~ 0(10^) at high temperature instead of 0(10^). The discrepancy is due 
to rates and asymmetries that occur in the Bohzmann equations and that are shghtly difTerent from the usual rates 
and CP-asymmetries in section [3] due to the different statistical factors they employ, for example 1 in equation (|79p 
or (1 — 2/]^') in equation (|92|) instead of the usual factor (1 — /^) in equation (|23p . 
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Figure 12: Evolution of neutrino abundance Yi^{z) and lepton asymmetry Yc{z) for K = 0.005 and 
zero initial neutrino abundance. We show the the two-mode cases and the vacuum case. 
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Figure 13: Evolution of neutrino density and lepton asymmetry for K 
neutrino abundance. We show the one-mode cases and the vacuum case. 



0.005 and zero initial 
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lepton asymmetry in the minus-mode obeys the same power law as the other finite-temperature 
modes, but is about a factor 100 lower due to the lower rates. The combined (ib)-abundance closely 
follows the plus-abundance since the influence of the minus-mode rates is very suppressed and also 
the different washout term with factor 1/4 instead of 1/2 can be neglected. 

Before turning to the intermediate temperatures z ~ Zgqj let us discuss the low temperature 
regime. For z > z^q, Iat » so the source term dominates and washout can be neglected. Since 
in this regime, Df^^ ~ D^^^ ~ D^^^^ we can write 



•^^^ r,(v v^^\ 
— ^ eoDiY^ - V) = -e— 



e,D{Yr, - F^^) = -6^ . (114) 



To first order, the negative lepton asymmetry created below Zeq and the positive contribution from 
above Zeq have the same magnitude and cancel each other. For the remaining asymmetry that did 
not cancel, the washout contribution up to Zeq and the exact behaviour of the abundances around 
Zeq are crucial. 

Assuming that Yp^[z = oo) = 0, we get 

y£fi"~eoY/v(^eq)-|l£(^eq)|, (115) 

so we see that the evolution of the difference AY{z) = Y]\f{z) — \Yc{z)\/eo below Zeq is crucial for 
the final lepton asymmetry. For the regime 1 ^ ^ z^q, we can write 

z 1 / Aj\ 



For the finite-temperature cases, the CP-asymmetry in the decay rates, A^/eq, is considerably 
smaller than the decay rate in the range Zq < z < 2, which can be seen in figure [HI Moreover, 
for the one-mode cases, A^ym/^o approaches faster than approaches 7+ for the plus-mode. 

Above z^, the ratio Aj/leoji)) for the two one-mode cases is about the same. Thus, the difference 
AD = D — is largest for the two-mode approach, smaller and about the same for the two 
one-mode approaches and vanishing for the vacuum approach, AD^ > AD^^ ~ ^/)v^"^f and 
AD^ = 0. As a result, AY^ > AY"^^ ~ ^yv^m^ therefore the final asymmetries are 

related as Y^ > Y^^ ~ y^™^ ^ Y^. Note that the final asymmetry is non- vanishing for the 
vacuum case, since the washout at higher temperature is larger than at lower temperature due to 
the larger decay rate. For the finite temperature cases, the difference D — in the crucial regime 
z ~ Zcq governs the final asymmetry. 

The evolution of the decoupled minus-mode at low temperature is very different but not hard 

to understand. The asymmetry rises at 2: > z^, because this is the regime above Zc where A7 ^. 

is maximal and therefore D~ is maximal as well. Between the thresholds and z^ , A7 (_ 

is suppressed by the internal plus-lepton and A7 is suppressed by the residue of the internal 

minus-lepton, so Z)^ is negligible and Y^ does not rise. Above 2; ~ 1, A7 ^. falls due to the residue 

of the external minus-mode and Y^ does not change. The final lepton asymmetry therefore does 
not change its sign above z > Zcq and keeps the value it achieves at around 1 ^ 2; ^ 2 where the 
CP-asymmetry A7 |_ becomes small. 

The combined (ib)-mode does not evolve differently from the plus-mode since the influence 
from the 7_ rates can be neglected and also the washout term with the additional factor 1/2 is 
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Figure 14: Evolution of neutrino density and lepton asymmetry for K = 0.005 and zero initial 
neutrino abundance. We display the four modes from figures [12] and [13] that give different final 
lepton asymmetries, that is, the plus-mode the minus-mode, the m^-mode and the vacuum case. 

not noticeable since washout is very small in all temperature regimes. Summarising, there are 
four differing lepton asymmetries in this regime: the vacuum case, the m^-case, the (+)-case and 
the (— )-case. The \/2m^-case yields the same asymmetry as the m^-case and the (±)-case yields 
the same asymmetry as the (+)-case. We show the four differing lepton asymmetries together in 
figure [H] 

Strong and intermediate washout for zero initial abundance 

For strong washout, the evolution of the neutrino abundance is analogous to the weak washout 
regime, with ~ z and ~ z^, as shown in figures [15] and [161 The couplings are stronger, 
therefore the abundances rise faster and meet Y^ earlier at Zeq ~ 1- For larger z, the couplings are 
strong enough to keep Yat close to equilibrium. The evolution of the lepton asymmetry is nicely 
explained in reference |49j for the vacuum case with some rather accurate analytical approximations. 
In this work, we are only interested in the difference of the vacuum case to the finite temperature 
case. In the strong washout regime, the lepton asymmetries rise rather fast and the washout term, 
which competes with the source term, becomes larger than the latter at some temperature z^i^, 
where the lepton asymmetry reaches its most negative value. The source term becomes small when 
Yn approaches its equilibrium value, so the washout term drives the asymmetry evolution back to 
zero. At z > Zeq, the neutrino abundance slightly overshoots Y^, so the source term changes sign 
and adds to the washout term until the lepton asymmetry becomes positive, where the washout 
term changes its sign as well and is competing again. At low temperature, source term and washout 
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Figure 15: Evolution of neutrino abundance Yi\i{z) and lepton asymmetry Yc{z) for K = 100 and 
zero initial neutrino abundance. We show the two-mode cases and the vacuum case. 
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Figure 16: Evolution of neutrino abundance Yi\f(z) and lepton asymmetry Yc{z) for K = 100 and 
zero initial neutrino abundance. We show the one-mode cases and the vacuum case. 



29 



term have the same magnitude when the lepton asymmetry reaches a maximum at Zmax- Above 
-Zmax) the lepton asymmetry is again driven to zero by the larger washout. At very low temperature, 
washout and source term become very small and do not influence the asymmetry further, which 
settles at a final value Y^^. We see that in the strong washout regime, the dynamics are governed 
by the washout term. The evolution of the finite temperature lepton asymmetries is analogous 
to the vacuum case, but they settle to a different final value. For the finite temperature cases, 
the equilibrium density of leptons is smaller than in the vacuum case due to the thermal mass 
nil ~ 0.2 T. The washout term is effectively larger than for the vacuum case and competes with 
the source term in a stronger way. Therefore, the asymmetry evolution appears slightly damped 
compared to the vacuum case and the final asymmetry is marginally lower. The evolution of the 
minus-mode is analogous to the evolution at weak washout, rises fast below and does not change 
above the thresholds since Yat ~ in this regime. The combined (ib)-mode tracks the plus-mode 
until washout becomes relevant at z > Zeq. The washout term for the (ib)-mode is always about a 
factor two smaller than for the plus-mode, since we add the minus- and plus-washout rates, where 
the minus-rate is always negligible compared to the plus-rate. Thus, the (it) -abundance is less 
affected by washout, so the dynamics are affected by the source term in a stronger way and the 
final asymmetry is larger. We can view this behaviour as always distributing half the asymmetry 
in a mode i- which couples strongly to and is not affected by washout. The final asymmetry is 
about a factor two larger than for the other scenarios in the strong washout regime. 

The case of intermediate washout is shown in figure [TTl where we only show the two- mode cases 
since in this regime, the final lepton asymmetries of the one-mode cases are the same as for the 
plus-mode. The dynamics can be viewed as an interpolation between the strong and weak washout 
regimes and the final asymmetries are very similar to each other. 

Non-zero initial abundance 

We also present the dynamics for thermal and dominant initial abundance in figures [T8VI23I For 
weak washout and thermal initial abundance, Yat ^ \Yc\ for low temperatures, and according 
to equation (IllSp . Y^'^'^'/^o ~ Y^'^''^^ For weak washout and dominant initial abundance, this 
equation holds as well, as can be seen in figures [18] and [T9j For intermediate washout ~ 1 and 
dominant abundance, shown in figures [20] and [211 the lepton asymmetry production is stopped 
between the thresholds for the thermal cases and the production above z ~ 1 does not succed 
in producing an asymmetry as high as in the vacuum case. For the (ib)-case, the asymmetry 
production is larger since it is not as much affected by washout. 

For the strong washout regime and large initial neutrino abundances, the dynamics at high 
temperature are interesting, as shown in figures [22] and [23l but the interplay between source term 
and washout term at low temperature governs the final asymmetry as in the zero-abundance case. 
We reproduce the well-known fact that the initial conditions do not infiuence the final asymmetry in 
the strong washout regime, while the arguments concerning the equilibrium distribution of leptons 
with thermal mass and the reduced washout of the (ib)-mode still hold and lead to the same lepton 
asymmetry as for zero initial abundance. The decoupled minus-mode is very much affected by the 
coupling, that is the decay parameter K, and the initial conditions, since the final lepton asymmetry 
is produced at high temperatures. The stronger the coupling, the larger the asymmetry production 
of the minus-mode at high temperatures and the larger the final value. Moreover, the larger the 
initial deviation of the neutrino abundance from equilibrium, the larger the asymmetry production 
and the final asymmetry. The final lepton asymmetry in this mode is thus lowest for neutrinos 
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Figure 17: Evolution of neutrino abundance Yi\i{z) and lepton asymmetry Yc{z) for K = 1 and 
zero initial neutrino abundance. We show the two-mode cases and the vacuum case. 
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Figure 18: Evolution of neutrino abundance Yn{z) and lepton asymmetry Yc{z) for K = 0.005 and 
thermal initial neutrino abundance. We show the two-mode cases and the vacuum case. 
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Figure 19: Evolution of neutrino abundance Y^[z) and lepton asymmetry Yc{,z) for K = 0.005 and 
dominant initial neutrino abundance. We show the two-mode cases and the vacuum case. 
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Figure 21: Evolution of neutrino abundance Yis[{z) and lepton asymmetry Yc{z) for K = 1 and 
dominant initial neutrino abundance. We show the one-mode cases and the vacuum case. 



with thermal initial abundance. 
Final lepton asymmetries 

The values of the final asymmetries are shown in figures [24H26] for different initial abundances. 
For zero initial abundance and weak washout, shown in figure [Ml the asymmetries for the finite- 
temperature cases are larger than for the vacuum case by about one order of magnitude due to the 
difference of the thermal rates 7£) and the CP- asymmetries A7 at z > . The lepton asymmetry 
for the plus-mode is also slightly larger than for the one-mode cases due to a suppression of the CP 
asymmetry compared to the one- mode approaches. For strong washout, the asymmetry production 
in the vacuum case is marginally more efficient than in the thermal cases due to a smaller lepton 
equilibrium distribution, while the lepton asymmetry in the (zb)-approach is by a factor two larger 
than in the other cases since half of the asymmetry is stored in the ^_ -modes and not affected by 
washout. The minus-mode case is completely decoupled, the lepton asymmetry bears the opposite 
sign as the other lepton asymmetries and rises with stronger couplings, that is with larger decay 
parameter K. As discussed in section l4?3l this scenario might not be realistic since the modes will 
couple to each other via gauge bosons, so an evolution similar to the (ib)-case seems more likely. 

For thermal initial abundance and weak washout, shown in figure [25l the final asymmetry 
equals the equilibrium abundance while in the strong washout regime, it shows the same 

behaviour as in the case of zero initial neutrino abundance. The minus-mode asymmetry is very 
low for thermal initial neutrino abundance since the neutrinos are close to equilibrium at high 
temperatures. Contrary to the zero initial abundance case, it bears the same sign as the lepton 
asymmetries of the other scenarios. 



33 



Yx 10-5 




^ N 



0.06 0.1 z. 



Figure 22: Evolution of neutrino abundance Yi\i{z) and lepton asymmetry Yc{z) for K = 100 and 
thermal initial neutrino abundance. We show the two-mode cases and the vacuum case. 
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Figure 23: Evolution of neutrino abundance Yn{z) and lepton asymmetry Yc{z) for K = 100 and 
dominant initial neutrino abundance. We show the two-mode cases and the vacuum case. 
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Figure 24: Final value of the lepton asymmetry for different values of K for zero initial neutrino 
abundance. 
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Figure 25: Final value of the lepton asymmetry for different values of K for thermal initial neutrino 
abundance. 
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Dominant initial abundance 
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Figure 26: Final value of the lepton asymmetry for different values of K for thermal initial neutrino 
abundance. 

For dominant initial abundance, shown in figure [26l the final lepton asymmetries assume their 
maximal value in the weak washout regime, when the coupling is weak enough not to wash them 
out at low temperature. For larger couplings ~ 1, the thresholds lead to a halted asymmetry 
production for the finite temperature cases and not as much asymmetry can be produced as for the 
vacuum case. The (ib)-case shows a larger asymmetry compared to the other thermal cases due to 
the weaker washout. At strong coupling K ^ 1, the asymmetries are the same as for thermal and 
zero initial neutrino abundance. The minus-mode asymmetry is large in all washout regimes, since 
the neutrinos are far from equilibrium at high temperatures when the £_-asymmetry is produced. 

A summary of the several initial conditions can be seen in figure [271 where we have omitted 
the '•/2 ml case since it is very close to the case in all scenarios. In the weak washout regime, 
the case with zero initial abundance is most strongly affected by thermal corrections which amount 
to one order of magnitude, and the plus-mode asymmetry is additionally enhanced by a factor 
of about two. In the intermediate regime, the dominant-initial-abundance case is influenced very 
much by thermal corrections. Therefore, a production mechanism for dominant neutrino abundance 
has to take into account such thermal effects. In the strong washout regime, one would naturally 
expect that thermal corrections can be neglected. We see that this is not the case since for strongly 
interacting leptonic quasiparticles, a part of the lepton asymmetry can be hidden in the ^_-mode 
which is unaffected by washout, thus producing an asymmetry by up to a factor of two larger than 
at zero temperature. The effect of the thermal lepton mass on the equilibrium distribution of the 
leptons is an interesting feature, but very small and can be neglected for all practical purposes. 
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Zero, thermal and dominant initial abundance 
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Figure 27: Final value of the lepton asymmetry for different values of K for zero, thermal and 
dominant initial neutrino abundance. The letters z, t and d denote the curves for zero, thermal 
and dominant abundance. Note that the final asymmetry of the minus-mode has opposite sign for 
zero initial neutrino abundance, compared to the asymmetries of all other cases. 
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5 Conclusions 



For a minimal and self-consistent toy model of leptogenesis, which consists only of neutrinos, leptons 
and Higgs bosons, we have performed an extensive analysis of the effects of HTL corrections. This 
implies capturing the effects of thermal masses, modified dispersion relations and modified helicity 
structures. We put special emphasis on the influence of the two fermionic quasiparticles, which 
show a different behaviour than particles in vacuum, notably through their dispersion relations, 
but also the helicity structure of their interactions. Our work is thus similar to the work done 
in reference [15], where the authors of the latter work did not include the effects of fermionic 
quasiparticles and get a different result for the CP- asymmetries, which are crucial for the evolution 
of the lepton asymmetry. Our toy model produces two lepton asymmetries stored in the two 
different lepton modes without the possibility of an equilibration of these asymmetries by SM 
processes. Since we expect the lepton modes to interact via gauge bosons in the bath, we examine 
a second case where the modes are strongly coupled to each other. As a third and fourth case, we 
approximate the lepton propagators by zero temperature propagators with the zero temperature 
mass replaced by the thermal lepton mass in one case and the asymptotic mass in the other case. 
We refer to these cases as one-mode approach. All four thermal cases are compared to the zero- 
temperature case. 

We have calculated interaction rates and CP-asymmetries in references [21] , |25j and |23j , where 
a detailed analysis can be found. We present the rates and CP-asymmetries shortly in section [3j 
Neglecting the zero-temperature fermion mass, the resummation of HTL fermion self-energies re- 
sults in an effective fermion propagator that does not break chiral invariance and is split up into 
two helicity modes. The external fermion states therefore behave conceptually different from the 
ones with chirality-breaking thermal masses that have been inserted in the kinematics by hand. 
Moreover, one has to take care of one additional mode, which has implications for the Boltzmann 
equations. 

We derive and evaluate the Boltzmann equations in section performing the crucial subtraction 
of on-shell intermediate states in appendix [QPI . We compare the results of the Boltzmann equations 
for our five cases, that is, decoupled lepton modes, strongly coupled lepton modes, the one-mode 
approach with mi, the one-mode approach with and the vacuum case. We assume three 

different initial values for the abundance of neutrinos: zero, thermal and dominant abundance, 
motivated by different scenarios for the production of heavy neutrinos after inflation [15]. In the 
weak washout regime, we flnd that using thermal masses enhances the flnal lepton asymmetry by 
about one order of magnitude for zero initial neutrino abundance. This is due to the fact that the 
CP-asymmetry and the decay rate evolve differently at z > 1 when using thermal masses, since 
the CP-asymmetry suffers from an additional suppression by thermal masses through the leptons 
and Higgs bosons in the loop. Due to the helicity structure of the modes, the CP-asymmetry 
of the plus-mode is additionally suppressed, which results in an additional enhancement of the 
plus-mode lepton asymmetry compared to the final asymmetries of the one- mode approaches. The 
enhancement we find is similar to the one found in reference [15j in this regime, but hard to compare 
quantitatively due to their different approach, which includes scatterings, and the discrepancy in 
the CP- asymmetry. 

In the strong washout regime, thermal masses do not show an influence, as expectecf^. However, 

^Reference [46] uses a thermal factor (1 — /jv) without explicitly deriving this factor. We show in this appendix 
that we have to use the equilibrium distribution for the neutrinos in (1 — /^) instead. 

^"There is a slight suppression of the lepton asymmetry for thermal masses, since the thermal mass suppresses the 
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when we couple the plus- and minus-mode strongly, we observe an enhancement of the lepton 
asymmetry by a factor of about two, since we stored half of the asymmetry in a mode that essentially 
does not interact with the neutrinos and is therefore not affected by washout. For intermediate 
washout, that is X ~ 1, we find that the lepton asymmetries with thermal masses are about one 
magnitude lower than in the vacuum case when we assume dominant initial neutrino abundance. 
This is due to the fact that the lower CP-asymmetry does not succeed in producing as much lepton 
asymmetry at 2; > 1 when using thermal masses. 

A decoupled minus-mode would show a behaviour completely different from the other thermal 
cases and the vacuum case for all initial values of the neutrino abundance. The lepton asymmetry 
in such a decoupled mode is produced mainly at high temperature and only slightly affected by the 
development at z > 1, where it decouples from the evolution of the other abundances. Therefore, 
the washout parameter K, which determines the coupling strength and thereby the asymmetry 
production at high temperatures, is crucial for the final value of the lepton asymmetry in this 
mode, as is the initial neutrino abundance. 

Summarising, we argue that for an accurate description of medium effects on leptogenesis, 
the influence of thermal quasiparticles, notably the effects of the two fermionic modes, cannot be 
neglected. Similar to reference our study shows that thermal masses have a strong effect in the 
weak washout regime, while the effect of fermionic modes has an additional influence on the final 
lepton asymmetry in this regime. We also showed that notably in the strong washout regime, the 
presence of a quasi-sterile lepton mode that is not affected by washout can have a non-negligible 
effect on the final lepton asymmetry. Future studies should clarify the dynamics of the interaction 
between the two fermionic modes and determine whether the evolution of the asymmetries in the 
two modes is closer to the decoupled or the strongly coupled case. 

Another important aspect that might be studied in future works is the influence of the finite 
width of the fermionic modes [51], notably the minus-mode. Such effects could be studied using 
formalism that takes into account non-equilibrium quantum effects, such as Kadanoff-Baym equa- 
tions p ^[T8H20| l22 |[52H57] . In the quest for a unified description of finite-temperature effects on 
leptogenesis, it is important to include SM interactions in the Kadanoff-Baym studies that are 
under way. To this end, quasiparticle excitations of fermions and gauge-bosons should be taken 
into account. Last but not least, the fermionic modes might have an influence on other related 
dynamics in the early universe that involve fermions, such as thermal production of axions, axinos 
or gravitinos, which could be studied in future works. 
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A Particle Kinematics 

The Boltzmann equations describe the time evolution of the distribution function of a particle 
species ^. We assume an isotropic and spatially homogeneous universe described by the Friedmann- 

equilibrium distribution of the leptons somewhat and thereby enhances the washout term. 
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Lemaitre- Robertson- Walker (FLRW) metric [47] . 

ds^ = dt"^ - a{tf I + ^^^^^ + ^^"^^ ' (^1"^) 

where a{t) is the cosmic scale factor, which describes the expansion of the universe, k = ±1,0 
specifies the curvature, and {t,r,9,(p) are the comoving coordinates. 

The trajectory of a particle ip with mass > moving in a gravitational field is given by the 
geodesic equations of motion [58] : 

^ + r!:„p>5 = o, (118) 

Since s = m^r is the eigen-time of the particle, r is fixed and is the momentum of a particle ■0- 
In the FLRW metric the /i = component of Eq (jllSp is given as 

Writing p° = |p^| d|p^|, this leads to: 



|p^|a + a|p^| = 

^ Pi/i = const. X -. (121) 
a 

Therefore the 3- momentum scales as 1/a. 

In general, the Liouville operator describing the evolution of a point particle's phase space in a 
gravitational field is given by 

L=p''^ T$yp^—. (122) 

With this operator the equations of motion (jllSp and (jllOp can be written for the momentum as 

^ = L [p^] , (123) 



and for the space-time as 



L [x^"] . (124) 



dr 

Furthermore, the time derivative of the phase space distribution of a non-interacting gas vanishes, 

= 0. (125) 

dr 
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Using the equations of motion for the particle we obtain the Boltzmann equations for the non- 
interacting particle species ip, 



L[U{x,p)]=0. 



(126) 



Since we assume a Robertson-Walker universe which is isotropic and homogeneous, the distribution 
function depends only on t and |p^|. Therefore, the Boltzmann equation can be written as [47] 



' dt 



! 



0, 



(127) 



where we have omitted arguments for the sake of notational clarity. 

Since = and because of the spatial isotropy of the Robertson- Walker-Metric, we have 



-E'V' I Pi/i I 



After dividing by E^, equation (|127p has the form 



L' [U] 



dt 



I ^■^V' 
d\P4>\ 



(128) 



(129) 



Interactions are introduced on the right-hand side by a collision term C[/^], which drives the 
distribution function towards its equilibrium value. The complete Boltzmann equation reads 



dt 



c[U]. 



(130) 



Thus, the Boltzmann equation in a Robertson- Walker universe has the form of a partial differential 
equation. However, in the radiation dominated phase of the universe, in which leptogenesis takes 
place, equation (I130p can be written as an ordinary differential equation by transforming to the 
dimensionless coordinates z = m^/T and = \p^\/T. Using the relation dT/dt = —HT, the 
differential operator dt — \p^\Hd\p.\ is written as zHdz, and consequently [59] 



dz 



fisil 



with H (m^) = H |T=m^ • In this form, the Boltzmann equation can be solved numerically on a 
grid for specific rescaled momenta y. For the right hand side, we have to sum over the collision 
terms of all processes which involve the particle '0 ^-nd change the ph ase space distribution. The 
collision term for a process ^/J + a + ■■■■<r^i + j + ■■■ is given by 



C[ii + a + 



+ = j W_dpa{'iT^f5^{p-^+Pa + Pi 



Pi 



) 



\M{i^ + a + ■ ■ ■ ^ i + 3 + ■ ■ - r Ufa- ■ ■ {l± fi){l± f,) 
|A1(i + j + • • • ^ + a + • • • )l Vi/i • • • (1 ± ± /a) • 



(132) 



^We have chosen a normahsation different from Kolb and Turner, so Chcrc = o-^Ckt 
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where a = (a, • • • , i, j, • • • ), 

= (2^)lk- ^'''^ 

The terms (1 ±/j) hold for fermions (— ) and bosons (+) and are interpreted as Fermi-blocking (— ) 
and Bose-enhancement (+). In practice, we will only look at processes which involve three or four 
particles, that is, decays, inverse decays and scatterings. We have included the internal degrees of 
freedom, g^,ga,- " ,9i,gjr"^ i^i the matrix elements, therefore we need to put in front of the 
collision term since it is not included in the phase-space density 

We integrate equation (jl3ip over the phase space of the incoming particle with g^ J d^p^/(27r)'^ 
and arrive at 



dn. 

7^ 

; Hirrinh) 

processes 

where 



H7 = "m^ ^ [7(^^ + + •••^i + j + ---)-7(^ + j + ---^V' + a + ---)], (134) 



7(^ + a + yi+ j + ...) = -g^ I l^C7[^ + a + + j + 



= / n dP/3(2vr)'^5'^(j>^ +Pa-\ Pi-Pj ) 

X \Mii^ + a + ■ ■ ■ ^ i + j + ■ ■ ■)\^ Ufa - ■ ■ {1± mi± fj) 



(135) 



where we now integrate over as well, that is, /3 = {ip, a, • • • The analogous equation 

holds for 7(i + j + -- -^'ip + a + -- -). 

B Boltzmann Equations at Zero Temperature 

We can derive the Boltzmann equations for the neutrino and lepton evolution at zero temperature, 
approximating the phase-space densities with Maxwell-Boltzmann distributions, 

fiiEi) = expi-Ei/3) , (136) 

where energy conservation in scatterings and decays implies 

fN = hfH (137) 

and there are no Higgs decays at high temperature. For the neutrino evolution, we get, analogous 
to equation (|i9]) . 

^ = -^(x7v-l)7o, (138) 

where 



70 



j dpNdpLdpHi27r)''6\pN - PH - PL) lA^ol' Q- (139) 
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The matrix element evaluated at zero temperature reads 

\Mof = 4.x2Pn-Pl, (140) 
where the factor 4 originates from summing over £ and i, as well as over the doublets (e~, (j)'^) and 

We can express the decay rate 70 in terms of the total decay width rj;^* in the rest-frame of the 
neutrino, 

where = 2 accounts for the internal degrees of freedom of the neutrino, the two spins, and 

Tlf{N^HL)=^-^^^^hlMl (142) 

describes the decay of a neutrino with a definite spin into (0^) and 
Evaluating equation (jl4ip . we get 

7o = gN^TTlfK,{z), (143) 

where z = M/T and Ki{z) is a Bessel function of second kind. For the equilibrium density of the 
neutrinos, we get 



so that 



70 ^-ptot -^i(^) 



(145) 



We can write the Boltzmann equation as 

Y^ = -D{Yn-Y;,'^), (146) 

where 



z 70 ^ ^i^jMf) 



D = —^ = zK^^, (148) 



^^'^ = ^--^^2(.) (149) 
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and 

-ptot ~ 

K = tj^ = ^ (150) 

is called decay parameter. 

For the lepton evolution, the subtraction of on-shell propagators can be performed analogously 
to the finite temperature case in appendix (Cj so that 

^sub(^^ ^ _ ^sub(|^ ^ ^ ^^^^^ (^5^) 

where 

is the CP-asymmetry. Analogous to equation ([67|) . we get 



which can be rewritten as 



where 



Y^ = eoD{YN-Y^'')-WYc, (154) 



We have 



where gt = 2 accounts for the lepton doublet components, so we get 

1 9iv .3 , 



W = -—z^KK2(z). (157) 
4 5£ 



C Subtraction of On-Shell Propagators 
C.l Low Temperature 



We verify the relation in equation (j62p . The scattering rate ^{ifp ~^ '^4') can be split up into four 
scatterings with different kinematics, corresponding to the four possibilities of combining the in- 
and outgoing lepton modes. The scattering rates read 



X \M{£H,<P^h^^)\''fiHMl-fgf,^)il + f^), (158) 
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where {hi,hf) = ±1 denote the hehcity-to-chh'ahty ratio of the initial- and final-state leptons (or 
antileptons). We will drop the subscript for this appendix part, unless it is necessary, and all 
equations are valid for one specific mode for each involved lepton, unless otherwise noted. With 
this simplified notation, each of the four matrix elements is evaluated as 



■ t 1 2 r 

{Vl)u \Dn\ 2 2{pN • PihXPN ■ Pihf) - {PN • PN){pth, ■ Plhf) 



(159) 



where we sum over the lepton spins S£ and and the lepton flavours and = '^/[P'n ~ ^-^n + 
ip^rAr(p^)] is the neutrino propagator in the narrow-width approximation and Tj^{p^) the total 
width of the neutrino, which equals the total interaction rate, including both lepton modes. Putting 
the propagator on its mass shell, = Mfj, we get 

2 
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2(pAr • PihdiPN ■ Pihf) - M^iPehi ■ Pihf) > (160) 



where 



7r<5(p2 _ m2) 



N\ 



P%^n{p%) 



(161) 



In vacuum without thermal masses, this reads 



I^^T2 



1 + 7?) 



(162) 



where the dependence on the angle rj between the external leptons cancels out in the integration 
for symmetry reasons, so we can neglect it and write 



\M{£c^ ^ A^)|^ |L>^f \MiN 



(163) 



At finite temperature with quasiparticle dispersion relations, we can not derive equation (jl63p 
accurately, but in the narrow-width approximation [15], one assumes that the influence of the 
angle between the external particles is negligible and equation ()163p holds. 
Using the relations in equation (j58|) . we derive 



\M°'{hA ^ 4,0/)! fih,U4^ - /ft,)(l + 4/) - \M°'{ihA ^ 4,0/)r /^-,,4i(i - fm,){l + Uj) 



\D 



N\ 



\M 



|2 



' [fcHM - + fZfcH, - 4e^/;,Vl - Q,)] Q{1 - fZ), (164) 



where we have neglected terms of order and xj^ and added the subscripts i and / in the Higgs 
boson distributions to clarify which momentum to use. 



(165) 



and likewise for 
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For the tree-level, CP-conserving amplitude, we have 



\M 



tree / 



(166) 



For the full amplitude |A^al=2| , the on-shell part is also dominant. Since it is CP-conserving, we 
write 



\Mal- 



■■'2\hih 



in; 



l-^AL=2lL/t 



OS 1 2 

n\ 



Ml 



(167) 



and we get 



|A^tree(^^^^ ^ 4,0)|'/ftj0(l " /,-,,)(! + /^) 

- |A^*'-^<'(4> ^ 4,</')|' /,-,,/^(l - fm,){l + U) 



\M 



AL=2\h^hf 



(168) 



Subtracting equations (jl64|) and (|168|) . we derive 



7™'(4,0 



7^"'(4> 



(169) 



Using the relations 



feqi-eq _ fS^t-i _ f eq , feq\ 

and - Q) = (1 + /!^)(1 - 



(170) 
(171) 
(172) 



(1 - + f?) = (1 - - f2 + /r)> 

f eq i-eq 

which hold for = ivgh + i^^), it is straightforward to derive 

7^"^4,</' ^ 4,<^) - 7'"^4> ^ 4,0) = 7''^^4,'/' ^ 4>) - 7^"^4,<^ ^ 4», 

7°^(L/.,iI ^ L;,^./!) = jTciLh.H ^ L;,,F) (173) 

Inserting 1 = / d^pAr/(27r)^5^(pAr — Pih^ ~P(l>) into equation p73p . again using the first relation 
from equations (jl69p and the expression for the total neutrino width, 



^n{jPn) 



h;=±l 



^Pn 

we arrive at equation 

[7^""(4> ^ 4,'^) - f "^4> ^ 4,0) 

= E [7'"^4,0 ^ 40) - f "^4,0 ^ 4. 



peq\ 



(174) 



dp^dm.dp<^(2^)454(p;v -m. -P*)^;^' \Miff,if;H^ - /: 

^ e'Hje^iLh^H ^ iV). 



eq\ 



(175) 
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C.2 High Temperature 

For the u-channel resonance at high temperature when Higgs bosons decay into neutrinos and 
leptons while the neutrinos are stable, we can derive a relation similar to equation (j62p . The width 
in the on-shell neutrino propagator is then not the decay rate but an interaction rate which accounts 
for the processes where the neutrino interacts with the medium, that is, H — t- NL and NL — t- H. 
This width acts as a regulator of the li-channel resonance. 

In the narrow-width approximation, the on-shell amplitude reads 



(176) 



where the on-shell propagator is the same as in equation (I16ip . but the width r^r is given by the 
kinematically allowed processes, H — )• NL and NL — t- H. 
Using the relations in equation (j87p . we derive 



2 r 



fckAi - ni) + Qjch, + KfZ(^ - f2,) - /;?/), (177) 



Analogous to equation (I168p . we derive 



\M 



tree 



(4.</> ^ 4,<^)r - /,-,,)(! + /^) 
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so that 

7^"'(4> 
Using the relations 



f "b(4,</' ^ 4, 



<^e!,iLh,H^L^,^H). 



Jlh>J(j) — J N \J Ih ^ J 4> )^ 



and - /;,^)(1 - /:^) = (1 + /!'^)/,\%^ 



which hold for uj^ = uj£h + wat, it is straightforward to derive 

7^"^4.'/' ^ 4,'^) - 7'"^4> ^ 4,</') = 7'"^4,'/' ^ ihA 



(178) 
(179) 

(180) 
(181) 
(182) 



(183) 



Inserting 1 = / d^pAr/(27r)'^5'^' 



Pat) into equation (jl83p . again using the first relation 



from equations (jl80p and the expression for the total neutrino width at high temperature. 



^n{p%) = / '^PLhfPH{'2^T^f5'^{PH - PLhf -Pn) 



Ml 



\J Lhf "T J H 



(184) 
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we arrive at equation ([88]) , 
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